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Abstract 

The result of a recent measurement of the size of the proton [R. Pohl et aL, 
Nature 466, 213] performed on the base of the muonic hydrogen spectroscopy 
turned out to be significantly different, by five standard deviations, from the 
results derived from the atomic hydrogen spectroscopy. This large discrepancy 
could come from the calculations of the Lamb shift in atomic hydrogen and 
muonic hydrogen. Here we show that there is a gap in the standard bound- 
state QED that may be the source of the discrepancy. This gap originates 
in the fact that within the framework of this theory the QED corrections are 
described in terms of the respective Green functions. The character of the time 
evolution of a system which should manifest itself in the general definition of 
bound states as stationary states of the system cannot be described in terms 
of the Green functions. We present a consistent way of solving the bound- 
state problem in QED starting from the condition of stationarity of the bound 
states. Formulae for the energies and the vectors of the states of one-electron 
(muon) atoms derived in this way indicate that the standard bound-state QED 
does not obey the exact description of the atomic states and, as a result, the 
Lamb shift obtained in its framework should be supplemented by an additional 
"dynamical" energy shift. It is shown that in this shift natural nonlocality of 
the electromagnetic interaction that in describing the S matrix and the Green 
functions is hidden in the renormalization procedure manifest itself explicitly. 
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Introduction 

The development of the theory of quantum electrodynamics was stimulated 
by high-precision hydrogen spectroscopy performed by Lamb and Retherford in 
1947 that showed a small splitting between the 281/2 and 2Pi/2 states known 
as the Lamb shift. Despite the theory of QED is the most well-tested, accu- 
rate, and successful theory in physics its mathematical foundation is not secure 
because of the still unsolved problem of the ultraviolet (UV) divergences. The 
renormalization theory that in QED is used for removing these divergences ac- 
tually allows one to pass over this problem but not to solve them. Richard 
Feynman said in this connection; "I think that the renormalization theory is 
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simply a way to sweep the difficulties of the divergences of electrodynamics un- 
der the rug" As is well known, in QED the ultraviolet divergences can be 
removed from the S matrix and the Green functions, but cannot be removed 
from quantities characterizing the time evolution of quantum systems, since reg- 
ularization and renormalization of the scattering matrix leads to the situation 
in which divergent terms automatically appear in the Hamiltonian and hence in 
the Schrodinger equation From this point of view the results of the recent 
muonic hydrogen Lamb shift experiment Q at PSI that have led to the proton 
charge radius significantly different, by five standard deviations, from the that 
derived from hydrogen atom spectroscopy 0-0] value seem not surprising. As 
Jeff Flowers of Britain's National Physical Laboratory wrote in Ref. [8] , if these 
experimental discrepancy is confirmed then "high- accuracy work such as that 
by Pohl and colleagues, not the high-energy collisions of giant accelerators, may 
have seen beyond the standard model of particle physics" . 

Here we show that the discrepancy may originate from the fact that the 
theory of QED fails in describing the time evolution of a system. For this 
reason one is forced to describe the bound states of QED systems in terms of 
the S matrix or the Green functions. However, in general bound states should 
be determined as stationary states which evolve in time according to the law: 

|vl/„(<)) = e-^"*|*„(0)). (1) 

If in QED the Hamiltonian, i.e. the operator of the total energy, were well 
defined, then from this definition it should immediately follow that the vectors 
of the stationary states and their energies are eigenvectors of the Hamiltonian 
H l^'n) = En I'I'n)- Howcver, this is not the case. The time evolution of QED 
systems remains uncertain even after renormalization. Actually this demon- 
strates the limitation of the existing theory of QED. The problem could be 
passed over, provided the bound states of composite systems in QED could 
be described in terms of the renormalizable S matrix and the Green functions. 
However, the discrepancy between the proton radii deduced from atomic hydro- 
gen spectroscopy and muonic hydrogen spectroscopy shows that this may not 
be the case: the QED corrections that played an important role in the determi- 
nation of these radii had been calculated by using the standard methods of the 
bound state QED in which the energy levels of composite systems are deter- 
mined by the positions of poles of the respective Green functions. This makes 
it necessary to investigate the problem more precisely. Since locality has been 
argued to be the main cause of infinities in QED, it seems natural to resolve this 
problem by introducing a nonlocal form factor into the interaction Hamiltonian 
density. However, it turned out that such an introduction of a nonlocal form 
factor results in the loss of covariance. The origin of this is the fact that the 
Schrodinger equation is local in time, and the interaction Hamiltonian describes 
an instantaneous interaction. But in relativistic quantum theory, if we spread 
the interaction in space, then we should spread it in time as well. Thus, for 
the introduction of nonlocality in the theory to be intrinsically consistent, a 
way is needed of solving the evolution problem in the case when the dynamics 
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in a system is generated by a nonlocal-in-time interaction. In Ref. [9| it has 
been shown that actually the Schrodinger equation is not the only dynamical 
equation consistent with the current concept of quantum physics, and the most 
general dynamical equation (GDE) consistent with these principles has been 
derived. Being equivalent to the Schrodinger equation in the case of instanta- 
neous interaction, this generalized dynamical equation allows one to generalize 
the dynamics to the case where the dynamics is governed by a nonlocal-in-time 
interaction. 

The GDE is shown to provide a consistent way of solving the bound-state 
problem starting from the law ([T]) determining bound states. In this way for- 
mulae are derived that determine the energies and the vectors of the states of 
one-electron (muon) atoms. The values of the atomic energy levels determined 
by these formulae differ from the values determined by the positions of the poles 
of the Green function of the electron in the Coulomb field. This difference that 
will be referred to as the dynamical shift is the reflection of the fact that there 
is not the exact correspondence between the atomic states and the poles of the 
Green function. The dynamical shift is shown to get the contributions from the 
processes, in which the nonlocal nature of the electromagnetic interaction that in 
describing the S matrix and the Green function is hidden in the renormalization 
procedure manifests itself explicitly. 



1. Quantum mechanical rules and the generalized dynamical equation 

Paraphrasing Steven Weinberg the recent situation in QED may be char- 
acterized as follows. If it turned out some QED system could not be described 
by the theory of QED, it would be a sensation. However, if it turned out that 
the system did not obey the rules of quantum mechanics and relativity, it would 
be a cataclysm. Of course, the question is raised, in this connection, what are 
the basic quantum mechanical rules that must be satisfied in any theory. In 
the Feynman's book [ll| where a minimal set of physical principles that must 
be satisfied in any theory of fields and particles is analyzed, the only quantum 
mechanical principle included in this set is the principle of the superposition of 
probability amplitudes. This principle formulated as the result of analysis of 



the phenomenon of quantum interference 12| reads as follows. 

The probability of an event is the absolute square of a complex number 
called the probability amplitude. The probability amplitude of an event which 
can happen in several different ways is a sum of the probability amplitudes for 
each of these ways. 



In the Feynman formulation of quantum theory 1^ [13| this principle is used 
as a basic postulate. This postulate provides the general rule prescribing how 
to calculate probabilities in quantum theory, and can be used in different ways 
depending on the choice of the class of alternative ways in which events can 
happen. In the Feynman formalism the processes associated with completely 
specified paths of particles in space-time are used as such alternatives. The 
contribution from a single path is postulated to be an exponential whose (imag- 
inary) phase is the classical action (in units of h) to the path in question. Thus 
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in the case of such a choice of the class of alternatives the contribution of each 
of the alternative way must be specified from the every beginning. However, 
very surprisingly it has turned out Q that there is a much more general class 
of alternatives whose contributions need not to be postulated: it is enough to 
know a priori only the contributions from the alternative processes associated 
with a fundamental interaction while the contributions from other alternatives 
are determined by the requirement of the conservation of probability. This class 
of alternatives consists of the processes with completely specified instants of the 
beginning and end of interaction in a quantum system. With such a class of 
alternatives, the superposition principle allows one to represent the probability 
amplitude of finding a quantum system in the state [^2) at time t, if at time to 
it was in the state in the form [§| 



(*2|(7(t,to)|*l) = (*2|*l) + J dt2j dh{^2\S{t2,h)\'^l). (2) 

to to 

Here {^2\S{t2,ti)\^i) is the probability amplitude that the interaction in the 
system begins in the time interval (ti,ii + dti) and ends in the time interval 
(^2,^2 + dt2), and after the end of the interaction the system will be found in 
the state |^'2) (we use the interaction picture), if before the beginning of the 
interaction the system was in the state l^"!). By using the operator formalism, 
one can represent amplitudes < 5'2|C/(t, io)|^i > by the matrix elements of 
the unitary evolution operator U{t,to) in the interaction picture. The operator 
S(t2,ti) defined in the same way represents the contribution to the evolution 
operator from the process in which the interaction in the system begins at time 
tl and ends at time t2- As has been shown in Ref. for the evolution operator 
in the form Q to be unitary for any t and to the operator S{t2, ti) must satisfy 
the equation 

t2 t4 

{t2 - ti)S{t2,ti) ^ J dti J dtsiti - t3)Sit2,t4)Sit3,ti). (3) 

tl tl 

A remarkable feature of this equation is that it works as a recurrence relation 
and allows one to obtain the operators S{t2,ti) for any ti and t2, if S{t2,t[) 
corresponding to infinitesimal duration times t = t'2 — t[ of interaction are 
known. It is natural to assume that most of the contribution to the evolution 
operator in the limit t2 — >■ ti comes from the processes associated with the 
fundamental interaction in the system under study. Denoting this contribution 
by Hint{t2i tl) we can write 

S{t2,ti) H„,t{t2M)+0{T'), (4) 

t2— »ti 

where t = t2 — ti. The parameter e is determined by demanding that Hint{t2,ti) 
called the generalized interaction operator must be so close to the solution of 
Eq. ([3]) in the limit t2 —5- ti that this equation has a unique solution having the 



4 



behavior (g]) near the point t2 = ii- Actually, Eqs. ^ and ^ represent the 
quantum mechanical rules that must be obtained by any physical theory. Equa- 
tion ^ is the representation of the probability amplitude < 5'2|C/(i, to)|^'i > 
as a sum of the contributions from all alternative ways in which the event can 
happen for the chosen class of alternatives, and Eq. ^ is simply the unitarity 
condition in terms of these contributions. The wonderful feature of this rule is 
that it directly leads to the dynamical equation. If Hint{t2,ti) is specified, the 
"unitarity condition" allows one to determine S{t2,ti) for any ti and ^2 and 
hence to construct the evolution operator. Thus, being supplemented by the 
boundary condition the relationship ^ becomes the equation of motion for 
states of a quantum system. It is important that this equation is universal, and 
the specific features of a theory describing the dynamics of quantum systems 
manifest themselves only in the boundary condition Q and in the Hilbert space 
with which the theory is dealing. In the case when the quantum field theory 
is required for the description of physical processes the Hilbert space must be 
chosen in the form of the Fock space, and correspondingly the interaction op- 
erator Hint{t2,ti) should be constructed in terms of the field operators. In the 
case when the interaction operator is of the form 

i^^„^(^2,^l) = -2m5{t2-ti)Hi{h), (5) 

i.e., the fundamental interaction is instantaneous, the generalized dynamical 
equation Q turns out to be equivalent [9] to the Schrodinger equation with 
the interaction operator Hi{ti). A quantum field theory must also obey the 
relativity. This condition is satisfied in the case when the interaction Hamilto- 
nian density 'H{x) related to the interaction Hamiltonian Hj{ti) by the equation 
Hj{t) = J (PxHi{t = 0,x) is local, i.e. the interaction is local both in time and 
in space. But this locality leads to the UV divergences. At the same time the 
interaction operator of the form ([5]) is the specific case of operators allowed by 
Eq. ([4]). In general the interaction operators describe the interaction spread 
both in space and in time. In this case the dynamics of a quantum system is 
determined by the behaviour of the interaction operator in the Schrodinger pic- 
ture, H^j^j(a;) , in the limit of infinitesimal duration times r of the interaction (in 
the local case the dynamics is determined by the interaction operator specified 
at the point t — 0). Thus there is no finite scale of nonlocality in the case of 
such nonlocal interaction! In order to obey these rules the operator describing 
the time evolution in QED should be of the form © with S{t2,ti) satisfying 
Eq. ([3]), and hence the fundamental QED interaction should be described by 
one of the nonlocal interaction operators allowed by Eq. ^ and obeying the 
rule of relativity. Thus the problem of the UV divergences is reduced to the 
problem of finding the form of this operator. 

Formally Eq. ^ involves processes associated with vacuum-vacuum tran- 
sitions. Correspondingly the matrix element {'^2\U{t,0) of the evolu- 
tion operator defined by Eq. ([2]) is the product of the probability amplitude 
(^2! Uph{t, 0) l^*!) of the physical event under study and the probability ampli- 
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tude of the vacuum-vacuum transition: 



(VI/2I u{t, 0) 1*1) = (*2| Uph{t, 0) 1*1) (0| U{t, 0) |0) . 

This one-loop contribution is proportional to the space-volume V because of 
the translational invariance. For the same reason the contribution from the 
two-loop vacuum process is proportional to and so on. The amplitude 
(^'2! Uph{t,0) 1^1), of course, is a sum of contributions from the process with 
the specified moments of the beginning and end of the interaction 



(*2|;7ph(t,0)|*i) = (*2 1*1)+ J dt2 J dh (*2|5j,^(t2,il)|*l) 

to to 

However, in this case the processes of the interaction associated with the 
vacuum-vacuum transitions must not be taken into account. Correspondingly, 
Sph{t2,ti) must satisfy Eq. from the right-hand part of which one has 

to remove the terms proportional to V^, i.e., the terms associated with the 
processes which involve the vacuum- vacuum transitions: 



{h ~ ti)Sph(t2,ti) — / dti / dt^lti — t3)Sph{t2,t4:)Sph{t3,ti) — c.t. 



where c.t. stands for counter terms proportional to V^. 

In the further discussion we will use the general notation (5'2| U{t,0) 
and {"^2] S{t2,ti l^fi) for describing the "physical" processes and the GDE in 
the form ([3]) keeping in mind that the contribution from the above terms in this 
equation must not be taken into account. 

It is extremely important that, being the representation of the general quan- 
tum mechanical rules, Eqs. ([2]) and (jS]) allow one to obtain detailed information 
about physical processes without specifying the interaction operator. Such in- 
formation can be regarded as a direct consequence of the first principles. Let 
us now investigate, in this way, the bound-state QED. 

Expression ([2]) for (4'2| U{t, 0) in the Schrodinger picture can be rewrit- 
ten in the form 

00 

(*2|(7s(i,0)|^'i) = / d£;e-^^*(*2|G(£; + iO)|*i), (6) 
2vr J 

—00 

where 

G{z) = Goiz) + Goiz)T{z)Go{z) (7) 
with Go{z) = {z — Hq)^^ and 

00 

T{z) = i J d{t2 - h) exp[iz{t2 - h)] exp[~iHot2]S{t2,ti) exp[iHoh] (8) 
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In terms of the T operator defined by Eq. ([5]) the generalized dynamical equa- 
tion ([3]) can be rewritten in the form [§| 

dz ^ [z-En? 

where n stands for the entire set of discrete and continuous variables that char- 
acterize the system in full, and \n) are the eigenvectors of Hq. Correspondingly, 
the boundary condition @ takes the form 

<1^2|T(z)|Vl>, ,^ <V2|SWIV'l>, (10) 
\z\ — >■ oo 

POD 

B(z)=i dTexp{izT)H[:l{T), (11) 
Jo 

where 

H'lntir) = exp{~iHot2)HM{t2,h)exp{iHoh) (12) 

is the interaction operator in the Schrodinger picture. Equation ([6]) coincides 
in its form with the standard expression that relates the evolution operator to 
the Green operator. In the canonical formalism the Green operator is defined 
as the resolvent of the total Hamiltonian H: 

G{z) ^{z- H)-^ (13) 

and Eq. ([7]) is the definition of the T-matrix. On the other hand, Eq. ([6]) is 
simply one of the forms of Eq. ([2]) representing the superposition principle, where 
the Green operator G{z) which in general can be defined as the inverse Fourier 
transform of the evolution operator, given by Eq. ([6|), can be represented in the 
form ([7]). The Green operator defined in this way coincides with G{z) defined by 
Eq. (jl3p only in the case when the interaction in the system is instantaneous, 
i.e., in the case when the GDE Q is reduced to the Schrodinger equation. 
This definition of the Green operator is valid even in the case when the total 
Hamiltonian H makes no sense. This is very important because the Green 
operator G{z) characterizes the evolution of a system, and, in particular, the 
positions of its poles determine the energies of stationary states, and this makes 
room for solving the bound-state problem in QED in a consistent way. 



2. The bound-state QED and the Green operator 

In solving the bound state problem in QED it is convenient to include the 
Coulomb field into the free Hamiltonian from the every beginning. This leads us 
to the Furry picture in which the eigenstates \n) of the Dirac-Coulomb Hamilto- 
nian Hq [Hq \n) = En'^ \n)) are used as "free" states. The one-electron states 
'^n^^ corresponding to the discrete spectrum of energy are just the "bare" 
states of a one-electron atom, i.e. the atomic states, when the atom does not 
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interact with the vacuum. The real atomic states are results of dressing the 
bare states "^n"^ ^ by this interaction. As we show below, the GDE provides a 
new very effective way of solving this problem. 

In QED the operator T{z) describes not only the interaction between par- 
ticles but also their self-action. This problem can be overcome by reduction, 
which amounts to the propagator Go{z) describing the evolution of free particles 
being replaced by the propagator Go(z) describing the evolution of particles in- 
teracting with the vacuum and, accordingly T{z) being replaced by T{z), which 
describes particle interaction proper. Redefined in this way the free Green oper- 
ator Go{z) describes the evolution of the system in the case when the particles 
move freely or interact only with the vacuum. In this way the Green operator 
can be rewritten in the form 

G{z) = Go{z) + Goiz)M{z)Go{z) (14) 

From Eq. (|14p it follows that the "free" Green operator Go{z) should be of the 
form 

and the GDE can be rewritten in the terms of Cm{z) and the operator M{z): 
(m' I m) ^^^^ = - < m'\P^M{z) (Go(z))' A/(z)|to >, (16) 



AP^M{z) (Go(z)) M{z)\i 



d<^'\M(.)\nr> ^ _ ^ Tn'\F^,M[z) I Gn(z) I M(z)\m > 



- < m'|P„i^%^Go(z)Af(z)|m > - < m'|Of(z)Go(z)%^|m >, 

(17) 

and Pn + P,i^l. The 



where P„ is the projection operator on the state 
corresponding boundary conditions are as follows 

<m'\M{z)\m> <Tn'\P^B{z)\m>, (18) 

\z\ — >■ OO 

< m'\m > C,n{z) < m'\P^B{z)\m >, (19) 

2 I — > oo 

Now we can obtain the vector of the corresponding atomic state. The sta- 
tionarity of the bound state of an electron in the Coulomb field |5'„(i)) with the 
energy £"„ must manifest itself in the fact that the Green operator G(z) has a 
pole ai z = En and a residue being the projection operator on this state: 

G{z)= ^ +0(1), z^E^ (20) 

The positions of the poles of the Green operator of the form are determined 
by the positions of the poles of the operator Go (2) that in turn are determined 
by the equation 

^„-^i") -G„(S„) = 0. (21) 
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where n corresponds to the discrete spectrum. For the positions of these poles 
to determine the atomic states 1^'™), their residues should be the projections 
on these states, |^'„) (^'n|- In order to find these residues, one has to analyze 
all of the contributions to the Green operator having such a pole. Besides the 
operator Go{z), this singularity obviously takes place in 

P^Go{z)M{z) = F,Go{z)M{z)P^ 

and _ _ 

Af(z)Go(z)P„ EE P„^Go(z)Af(z)P„. 

Note that < *^°'|M(z)|* ^'i°^|M(z)F,^|* > and < ^\M{z)\^'^^ > = 

< ^\P^M{z)\^n^ > by the definition. The singularity in the part of the 
Green operator given by P.;^Gq{z)M{z)Go{z)P.^ can arise only because of the 
singularity of the operator P^M{z)P^. It is very important that Eqs. (|16p 
and (jl7p allow to determine the behavior of the operator T(z) in the vicinity 
of the singular point z = En without specifying the interaction operator, and 
this behavior is as follows 

Mi^) = [MiEn)\^^^^ X ^l^^\M{En) + P^M{z)Pn + PnM{z)P^] + 

+0(1), z-^En 

(22) 



where Co =(l-^|,^^^ 



-1/2 



The insertion of this expression in Eq. ([14)) yields 

I*' ) (* I 

Z £/„ 

with ^ 

\^n) = C, (|vl/(o)) + P^Gt{En)M+{En) (24) 

I*:.) = Co + P,iG,{En)M{En) (25) 

Pd z = El [n — 1 corresponds to the ground state of the atom) the operators 
Gq{z) and M{z) are Hermitian and, as a result, \^'^ = From this and 

the fact that the vector 1*1) is normalized | v&i) = 1, and Eq. ([231), it 
follows that the ground state \^i{t)) is stationary: |*i(t)) = e"'-^i* |^'i(0)). 
The difference SE^^ = C{Ei) between the energy Ei of this state and E\ is 
caused by the interaction of the atom with the vacuum, and is just the Lamb 
shift of the energy level of the ground state. The vector describes the 
ground state of the atom. This is not a one-electron state even in the sense of 
the Dirac equation. Actually in this case we deal with the state of the atom 
surrounded by a cloud of virtual particles. The second term in expression (j24l) 
for the state just describes the probability to find the virtual particles such 
as, for example, photons and electron-positron pairs in a measurement, if before 
the experiment the atomic system was in the state l^*!). 
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The exited states |'i'„(i)) (n = 2,3,) can be regarded as quasistationary 
states because their energies determined by Eq. ([2T|) are complex En = En^ + 



Cn{En) = E^n^ + SE^ - ?r„/2. This results in the fact that the states l^-^) and 
Y^n) given by Eqs. (|24p and (P5|l respectively do not coincide exactly. 

Our analysis of the unstable atomic states can be extended to the general 
description in terms of the energy distribution fl3 |. However, in this work 
we will restrict ourselves to the approximation that is proven to be applicable 
for solving the most of the problems in quasistationary atomic physics. In 
this approximation the difference between the vectors j^JJ and \^n) can be 
neglected. 

Formulae (|2ip and (|24|) determining the energies En and the vectors l^'n) f 
the atom have been derived as an inevitable consequence of Eqs. dTB]) and p?]) . 
which in turn represent the basic quantum mechanical rules, without resorting to 
any specific information about the character of the electromagnetic interaction. 
In order to use these formulae for the calculations one has to obtain, at first, 
the operator M{z) and the function C„(z) by solving Eqs. ([T6|) and ([T7| . Only 
at this stage of our analysis we have to specify the interaction operator. The 
interaction operator B[z) that describes the local electromagnetic interaction is 
of the form 

B{z) = Hj = J (fxHi{0,x) J d^x: ^^(0, x)*(0, x)7'^«'(0, x) : (26) 

where Hi{x) is the interaction Hamiltonian density, 5'(t,x) is the Dirac field 
in the Furry picture and A^{x) is the electromagnetic field. The value of the 
coupling constant of the electromagnetic interaction allows to solve Eqs. (|16p 
and (|17p perturbatively. In the first order of this solution we have 

M^'\z)^Hi, Gi'\z)^Go{z) (27) 

In the next order, for C„(z) we have the equation 

dCi?\z) 



dz 



<*W|i/z(Go(z))'ff/|vI/(")> (28) 



The form of the interaction operator (|26|) implies that C„(z) tends to zero as 
\z\ — oo. Solving Eq. ([28]) with this boundary condition yields 

oo oo _ 

{E) = / e^^"(*^-*i)C„(t2 - h)d{h -tl) = J d{t2 ~h)J d3pi / d3p2*„(p2) X 
X S^(p2, ^2; Pi, ii)e^^"(*^-*^) + S^(P2, ti; Pi, i2)eH(^-2^i°')(*^-*^)}l M.„(pi) 



(29) 



where 



SA(P2,i2;Pi,ti) = J d^xi J d^X2^A{x2,xi)e'P-^-e-'P'^' (30) 

with Y,a{x2, xi) being the ordinary one- loop self-energy operator of the electron. 
The leading order Cn{E) obtained in this way can then be used for obtaining 
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the leading order energy shift. For this one has to solve Eq. (|2ip which in this 
case is reduced to the equation 

E-E^^'^ -C^p{E + iO)^0 (31) 

Neglecting the dependence of Cn^ (z) on z in the vicinity of the point z = En'' , 
from this equation for the leading order self energy correction SEn = E — En'' , 
we get 

<5(°)£;„ = ci2)(£;i"^z0) = J d^pi J d3p2§„(p2)E^(p2,pi,i;(°))^'„(pi) (32) 

Here we have taken into account that I]a(P2, ^2; Pi, ^i) depends on the differ- 
ence (t2~ti) only. Of course, the one-loop self-energy operator T, a{P2, Pi, En'') 
makes no mathematical sense because of the UV divergences. But the renor- 
malization theory is applicable in this case, and the problem is solved by replac- 
ing T,AiP2, Pi, En^) in Eq. ([32]) by its renormalized value E^'' (p2, pi, -En°^). 
Eq. (|32p coincides with the well known expression for the one-loop energy shift 
in hydrogen that in the standard bound state QED is derived, for example, from 



the solution of the effective Dirac equation with the mass operator 15 1 



H^-i'ix) + i J -fo^A{x, x', E)^{x')dx' = E^{x) (33) 

where is the Dirac Hamiltonian without the mass operator but with the 
external field. In fact, in the leading order of the expansion of the solution 
of Eq. ([55]) in powers of a we can put 4'(x) — ^'0(2:), and T,a{x,x' , E) = 
E^(x, x.' , Eq) in this equation, and in this way we arrive at the expression p2p for 
the self-energy correction. However, ([5^ is not an exact one-loop energy shift, 
because E = E^^'' +Cn^ {En) is only an approximate solution of Eq. (PT|) . At the 
next order of the iterative solution of this equation we get z ~ En^ + C^^"* {E'n). 
By using Eq. ([29]) and the representation 



EA(x2,Xi;t2,ii) = 7^ f dEe''^^*^-''^j:A{^2,^i;E) (34) 
Ztt J 

at this order for the one-loop correction, we get 

SE^^^ = Jd^pi j d^p2^n{P2)^AiP2,Pl,E'J^n{Pl)+6^En, (35) 
S^'En^ J d^p, J d3p2^r„(p2)s(«)(p2,pi)^-„(pi) (36) 



where 



with 



e(")(p2,p0 = -|- / (37) 

2^*-^ (^E^^^-E + ^0^ 
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The first term on the right-hand part of Eq. psp is the energy shift, which is 
determined by Eq. ([32]) . when in this equation we put Sa(x, x', E) = Ea(x, x', Eq) 
Thus, there is a discrepancy between the predictions of the standard bound state 
QED and our approach. 

The discrepancy is the manifestation of the fact that actually there is no 
one-to-one correspondence between the poles of the Green operator and the 
poles of the respective Green functions: Eq. (34) determines the poles of the 
Green function of the electron in the Coulomb field. To explain this point note 
that the electronic Green function in the Furry picture is determined as 

G.fc(x,xo) =^(0|^*f)(x)*i^)(xo)5|0) (38) 

where S is the scattering matrix and "i['^\x) and "i^i^\xa) are the electronic 
field operators in the Furry picture, and i,j are spin indexes, in the one-loop 
approximation under study, the S matrix is given by the first two terms in the 
Feynman-Dyson expansion, and, as a result, the Green function takes the form 

G^k{x,xo) = G|^^(i,x;to,xo) = i (0| r*,(i, x)*^(to, xq) |0) - 

oo oo _ (nQ\ 

-| / / d^x2{o\T^s,^)'^k(.to,^a)ni{xi)niix2)\o) ^ ' 

— oo — oo 

On the other hand, the one-electron matrix element of the evolution operator 
< ^n^-'|C/(i,io)|^n"'' > hi this approximation reads 

t oo t oo 

{^„\U{t,to)\'^„)^l-'-Jdti J d'xi Jdt2 J d^X2 (*„|r?^/(xi)7^/(x2)|4'„) 

to —oo to —oo 

(40) 

Here as well as in Eq. (p9)) the contributions from the terms associated with the 
vacuum-vacuum transition should be removed. Thus, in order for the electronic 
Green function could be reduced to the one-electron matrix elements of the evo- 
lution operator, the integration in Eq. (|39|) must be limited to the time intervals 
to < ti < t and to < t2 < t. This proves that the positions of the poles of the 
electronic Green function do not coincide exactly with the positions of the poles 
of the Green operator, and hence do not determine exactly the atomic energy 
levels. Nevertheless, for atomic hydrogen the contribution from the virtual pro- 
cesses that should be removed by the above limitation is relatively small, and 
this is a reason why the standard bound-state QED provides high accuracy in 
calculating the Lamb shift of the atomic energy levels. But this accuracy may 
turn out to be not sufficient for the derivation of the proton radius. 

The above shows that actually there is a correction S^E to the radiative 
Lamb shift in hydrogen that is hidden in solving the problem in a standard way. 
But the integral in Eq. (|36|) that determines this correction in the one-loop 
approximation diverges at infinity, because the one-loop self-energy operator 
E^(p',p;£^) behaves in the limit \E\ — > cx) as 

. , X a , fm'^ - (E + iO)'^\ , , 

S^p',p;i? ^ ^7oi?ln ^ -] 41 

\z\ -J- oo 47r \ / 
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The problem is that this divergence occurs after renormahzation. 

It should be noted, the above results can be trivially extended to the de- 
scription of QED corrections to energy levels of muonic hydrogen. 

3. Nonlocality of the electromagnetic interaction 

We have shown that the standard methods based on the determination of 
the atomic levels by the positions of the poles of the Green function of the 
electron in the Coulomb field do not obey the exact description of the atomic 
states. The dynamical shift 5^ E which in the one- loop approximation is given 
by Eq. describes the part of the Lamb shift that was missing in the standard 
description and might be the origin of the discrepancy between the radii derived 
from atomic hydrogen and muonic hydrogen spectroscopy. At the same time, 
in describing this dynamical shift we face the problem of the UV divergences 
that cannot be cured by renormahzation. In other words, the renormahzation 
theory fails in this case. This is not surprising because in dealing with the 
Green operator constructed by using the renormahzation procedure we are go- 
ing beyond the domain of applicability of the standard theory of QED, which 
is limited to deal only with the process that can be described in terms of the 
S matrix or the Green functions. The "rug" has turned out to be too small to 
hide the problem of the UV divergences that arise in the consistent description 
of the bound states based on the employment of the Green operator. Actually, 
in this case the natural nonlocality of the electromagnetic interaction, which in 
describing processes associated with the S matrix is hidden in the renormahza- 
tion procedures, manifests itself explicitly. This means, that for the accurate 
calculations of the QED corrections and, hence, for the accurate derivation of 
the proton radius from spectroscopy of atomic hydrogen and muonic hydrogen, 
one has to deduce the nonlocal interaction operator describing the fundamental 
electromagnetic interaction. Keeping in mind that taking into account the non- 
locality of the electromagnetic interaction which must result in improving the 
behaviour of the self-energy operator at infinity, we can make some assumptions 
about the value of the correction S^E. Formally the order of magnitude of this 
correction is given by 6^En — a'^AoS'^^'' En, where S^^^En coincides with the 
ordinary one-loop Lamb shift. The constant Aq should be determined by the 
form of the self-energy operator. Thus the nonlocality of the electromagnetic 
interaction must manifest itself in this constant and, hence, in the Lamb shift 
of atomic energy levels. 

It might seem that there is a significant arbitrariness in choosing the form 
of the interaction operator. However, the fact that the interaction operator 
Hint{t2,ti) in the limit t2 ti should be close enough to the relevant solution 
of the GDE, imposes strong limitations on its form. In addition, the form of 
the interaction operator is constrained by the symmetries of QED and the fact 
that in describing the processes associated with the S matrix the GDE should 
lead to the same results that in the standard theory of QED are obtained by 
using the renormahzation theory. A transparent example of this feature of the 
GDE is given in Ref. [l6|, where it has been shown that after renormahzation 
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in the effective field theory of nuclear forces the low energy nucleon dynamics 
is governed by the GDE with a nonlocal-in-time interaction operator whose 
form is determined by the above constraints up to a constant that is fixed by 
fitting to the nucleon scattering data. Here we are not in position to discuss the 
problem of the determination of the form of the nonlocal interaction operator 
describing the fundamental electromagnetic interaction in detail. But note, 
that there is a hope the above limitations on the form of this operator will 
allow one to determine it up to a constant the value of which could be derived 
from the spectroscopy of atomic hydrogen and muonic hydrogen. And just the 
discrepancy between the radii derived from these measurements may provide 
us experimental information about the nonlocal nature of the electromagnetic 
interaction. 

4. Outlook 

We have shown that the GDE provides a consistent way of solving the bound- 
state problem in QED. In this way the energies and vectors of the atomic states 
are determined by the positions and residues of the Green operator defined by 
the inverse Fourier transform © of the evolution operator. From Eq. ([25]) it 
immediately follows that such states are stationary. Formulae (PT|) and (|24p 
derived in this way obey the accurate determination of the energies and the 
vectors of the atomic states. This is because they are derived as an inevitable 
consequence of Eqs. ([2]) and ([3]) that represent the basic quantum mechanical 
rules without entering into the details of the interaction in the system. As we 
have shown the standard bound-state QED does not obey the rigorous descrip- 
tion of the atomic states, because strictly speaking the problem of finding the 
poles and residues of the Green operator is not reduced to the problem of finding 
the poles of the Green function of the electron in the Coulomb field. 

The discrepancy between the proton radii deduced from atomic hydrogen 
spectroscopy and muonic hydrogen spectroscopy seems to be so serious that it 
could force revisions of the fundamentals of physics. However, a small part of 
what is allowed by the current fundamental principles of physics has been real- 
ized in existing theories. In fact these principles, as we have seen, allow nonlocal- 
in-time interaction operators as generators of quantum dynamics, while in the 
existing theories one restricts oneself to the instantaneous interaction operators 
which generate the Hamiltonian dynamics. And the UV divergences are just 
the cost for this restriction. This means that local interactions are actually in- 
consistent with the current physical principles. Thus in order to realize all the 
possibilities allowed by the current principles of physics, one has to take into 
account the nonlocality of the fundamental interactions explicitly, and the GDE 
provides a consistent way to succeed in this object. And what is important in 
this connection is that Eq. ([2]) and the GDE separated from boundary condi- 
tion Q are simply relations that represent the basic quantum mechanical rules. 
And these rules impose such strong limitations on the physical processes that 
Eq. ([3]) turns out to be the rule prescribing how to build the probability ampli- 
tudes of physical processes from the "bricks" associated with the contributions 
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from the processes with infinitesimal duration times of interaction. This rule is 
general and is independent of physics with which we deal. The specific physical 
content is introduced by choosing these "bricks", i.e. specifying the form of the 
fundamental interaction operator. In addition this interaction operator must 
obey the relativity and other symmetries of the theory. 

If the theory of QED were defined in this way, then it would be finite and 
could provide the accurate calculations of atomic energy levels and, hence, the 
accurate derivation of the proton radius from the hydrogen atom and muonic 
atom spectroscopy. However, additional experimental information is needed to 
specify the nonlocal interaction operator. And the work of Pohl and colleagues 
may provide this information, which being reanalyzed in this way may reveal 
unknown features of the electromagnetic interaction. And these features may 
turn out to be very surprising. 

I thank Aldo Antognini for the discussions which triggered this study. 
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